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Section  1 


riTROOUCTIOM 

This  report  is  the  final  documentation  oescribing  work  com- 
pleted under  Contract  OfJA  001-76-C-028^ . The  objective  of  this 
analysis  was  to  assess  the  rationale  for  additional  underground  tests 
(UGT)  to  produce  an  air  blast  environment.  To  accomolish  this  objec- 
tive, measures  relating  to  the  tactical  utility  of  achieving  Increased 
accuracy  in  predicting  nuclear  weapon  produced  environments  were 
developed  and  applied  to  existing  environments  data  (blast  and  radia- 
tion). This  work  was  accomplished  in  two  phases.  The  first  phase 
resulted  in  development  cf  the  methodology  to  assess  the  rationale 
for  additional  tactical  nuclear  underground  environments  testing,  and 
a preliminary  application  of  this  methodology  to  selected  weapon  pro- 
duced environment  data.  The  second  phase  was  oriented  toward  strength- 
ening and  broadening  the  theoretical  and  empirical  basis  for  the 
methodology,  testing  some  of  the  assumptions  upon  which  the  initial 
applications  were  based,  and  additional  applications  of  the  methodology 
for  a more  complete  assessment  of  the  UGT  rationale. 

This  report  summarizes  work  performed  in  phase  1,  and  presents 
a detailed  description  of  the  analysis  completed  in  phase  2.  The  work 
completed  in  phase  2 and  described  *'erein  includes  the  following: 

e Refinement  and  concise  synopsis  of  the  theoretical  and 
empirical  basis  for  the  rationale  assessment. 

• Testing  of  the  assumption  of  a lognormal  distribution  for 
the  weapon  produced  environments,  which  is  an  isSue  in  the 
application  of  the  methodology. 

• An  analysis  of  the  nuclear  blast  data  to  attempt  to 
ascertain  if  syste?nat1c  errors,  or  biases,  could  exist 
which  would  require  additional  measurements  of  weapon 
produced  environment  to  resolve. 


• An  analysis  of  the  peak  dynamic  pressure  blast  data  to 
assess  possible  testing  rationale  based  on  this  blast 
phenomenon. 
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In  assessing  the  rationale  for  continued  tactical  nuclear 
underground  environments  testing,  it  is  important  to  distinguish  in 
the  analysis  between  assuming  that  the  weapon  oroduccd  environment 
data  are  randomly  distributed,  and  assessing  the  rationale  based  on 
that  assumption  and  recognizing  the  possibility  of  the  existence  of 
systematic  errors,  or  biases,  in  the  environment  test  data  that  could 
have  ramifications  for  the  tactical  utility  of  the  weapons. 

In  an  operational  context,  the  effects  of  these  two  kinds  of 
errors  are  quite  different.  If  the  probability  of  kill  (P^)  is 
simplistically  viewed  as  a specified  weapon  produced  environment 
being  produced  at  a given  range  95  times  out  of  100,  then  this  notion 
can  be  easily  interpreted  in  terms  of  random  shot-to-shot  variation 
in  the  weapon  produced  environment.  Hoi'<ever,  if  a systematic  error 
causes  a bias  between  the  predicted  and  true  weapon  produced  environ- 
mcnts,  all  jyeracinnal  shots  could  be  -.ystematically  different  from 
predictions.  Thus,  instead  of  a specified  effect  occurring  at  a 
calculated  range  95  percent  of  the  time,  the  effect  could  never  occur 
at  the  predicted  range  in  an  operational  environment  if  there  was  an 
unmodeled  bias. 

Alternatively,  the  choice  of  yields  estimated  to  accomplish  a 
specified  objective  could  be  systematically  too  large  or  too  small 
because  of  such  a bias.  Thus,  this  assessmen  proceeded  along  two 
parallel  paths  with  regard  to  analysis  of  the  nuclear  blast  data, 
namely: 

• An  assessment  of  the  rationale  for  tactical  nuclear 
underground  effects  testing  based  on  the  assumption 
that  variations  observed  in  the  blast  data  are 
randomly  distributed. 

• A preliminary  analysis  of  the  existing  blast  data 
with  the  objective  of  attempting  to  identify  possible 
biases  (systematic  errors). 

The  blast  data  variations  would  be  classified  as  random  if  the  statis- 
tical nature  of  variations  observed  in  the  blast  test  data  were 
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assumed  to  be  the  same  as  the  variations  that  one  would  observe  during 
deployment  of  the  weapons.  On  the  other  hand,  biases  between  the 
statistical  characteristics  of  the  variations  observed  from  the  test 
data  and  during  weapon  employment  would  represent  a case  where  the 
predictive  models  are  either  under-specified  or  mis-specified. 
Obviously,  large  biases  could  adversely  impact  the  tactical  utility 
of  the  weapons. 

To  strengthen  the  random  error  analysis,  an  assessment  was 
made  to  ascertain  if  the  assumption  of  a lognormal  distribution  for 
the  weapon  produced  blast  environment  is  a reasonable  one.  The  random 
error  methodology  is  based  on  this  assumption.  This  approach  involved 
an  analysis  of  the  relationship  between  the  error  in  estimating  the 
range  to  a specified  nuclear  produced  environment,  and  the  range  for 
peak  overpressure  and  overpressure  impulse.  An  approximately  linear 
relationship  between  the  error  in  range  and  range  would  suggest  that 
the  assumption  of  lognormal  distribution  for  the  error  is  a reasonable 
one,  since  this  is  a characteristic  of  data  that  are  lognormal ly 
distributed.  This  analysis  and  the  results  are  described  in 
Section  4.2. 
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SUMMARY  OF  RESULTS  AMO  CONCLUSIONS 

The  results  of  this  uncertainty  analysis  indicate  that  the 
justification  for  improving  the  predictive  capability  of  the  nuclear 
produced  blast  environment  is  marginal  provided  the  variations  in  the 
existing  blast  data  from  nuclear  bursts  are  treated  as  random.  The 
quantitative  basis  for  this  result  is  that  the  target  vulnerability 
and  weapon  emplacement  uncertainties  essentially  dominate  the  weapon 
produced  environment  uncertainty.  (Conclusion  1} 

On  the  other  hand,  it  has  been  found  that  there  are  non- 
trivial systematic  errors  (biases)  in  some  blast  data  which  could  be 
included  in  the  predictive  models  of  the  blast  environment  if  the 
causal  physics  of  these  biases  were  better  understood.  The  impact  on 
this  uncertainty  of  one  or  more  additional  nuclear  tests  has  not  been 
quantitatively  evaluated  because  this  result  came  from  a secondary  task. 

The  general  methodology,  developed  in  Section  3 for  this  random 
uncertainty  analysis,  is  based  on  the  assumption  that  the  randomness  in 
the  nuclear  produced  environment  which  would  be  experienced  during 
employment  of  the  tactical  nuclear  weapons  is  statistically  the  same 
as  the  deviations  observed  among  past  weapons  tests.  The  methodology 
relates  the  magnitude  of  this  environmental  uncertainty  (potentially 
reducible  through  additional  underground  testing  of  nuclearly  produced 
environments)  to  a measure  of  the  tactical  utility  of  the  weapon  (i.e., 
yield  change).  Thus  the  direct  payoff  from  continued  underground 
detonations  of  tactical  nuclear  weapon^  in  order  to  more  accurately 
predict  the  environment, is  related  to  the  tactical  utility  gained  from 
this  improvement. 

To  obtain  closed-form  expressions  for  the  relationships  between 
the  measure  (yield)  and  the  environment  predictive  capability,  it  was 
assumed  that  the  uncertainty  in  the  nuclear  produced  environment  was 
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proportional  to  the  measurements  which  are  then  implied  to  be  distri- 
buted lognormally.  This  hypothesis  was  tested  in  Section  4.2,  and  the 
characteristics  of  the  uncertainty  in  the  blast  data  were  found  to 
support  the  assumption  of  lognormal ity.  (Conclusion  2) 

Using  the  above  mentioned  methodology  (which  assumes  that 
there  are  no  biases  between  the  weapon  produced  environments  observed 
during  testing  and  during  employment  of  the  weapons),  the  maximum 
possible  reduction  in  yield  (average  is  also  available  in  Section  4.1) 
is: 

• Peak  Overpressure:  21  percent 

• Peak  Dynamic  Pressure:  9 percent 

• Overpressure  Impulse:  15  percent 

It  is  concluded  that  only  a tenuous  rationale  exists  for 
underground  testing  of  tactical  nuclear  environments,  based  on  the 
assumptions  upon  which  the  rationale  was  assessed  (no  biases).  The 
dynamic  pressure  impulse  data  are  probably  too  sparse  and  too  uncertain 
to  be  statistically  analyzed.  Therefore,  in  cases  where  the  dynamic 
pressure  impulse  totally  dominates  the  target  damage  mechanisms,  there 
is  a rationale  for  additional  underground  testing  of  tactical  nuclear 
environments. 

The  analysis  !n  Section  5 for  the  existence  of  systematic 
errors  (biases)  was  conducted  fo*"  scaled  ranges  to  10,  30  and  ICO  psi 
peak  overpressure  data  from  tests  between  0 and  11  ft/kt^^^  scaled 
height  of  burst  (SUOB).  Such  biases,  if  found,  would  represent  blast 
predictive  modeling  under-specification  or  mis-specification,  and 
would  require  an  attempt  to  adjust  or  improve  the  predictive  models 
to  remove  the  bias.  Since  part  of  the  objective  of  tnis  bias  analysis 
was  to  assess  the  degree  to  which  the  Johnie  Boy  (negative  SII03)  test 
results  differ  from  surface  test  results,  the  analysis  was  limited  to 
surface  and  very  low  SHOD  tests. 
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The  conclusion  from  this  limited  bias  analysis  is  that,  for 

scaled  range  to  TOO  psi  and  30  psi  peak  overpressure,  there  are 

1/3 

statistically  significant  deviations  from  the  1/  scaling  law  over 
the  yield  range  between  15  t and  15  Mt,  a factor  of  one  million.  More 
precisely,  the  low  yield  data  shows  a deviation  from  cube-root  scaling 
of  high-yield  data  due  most  likely  to  a mass-to-yleld  effect.  This 
range  bias  is  relatively  small  but  not  insignificant  when  translated 
into  a potential  yield  mis-specification.  For  instance,  if  the 
medium  and  high  yield  data  were  used  to  estimate  scaled  range  to  30 
psi  for  small  yield  weapons  (or  vice  versa),  a yield  under-estimate 
(or  over-estimate)  of  approximately  25  percent  would  occur. 

(Conclusion  3) 

This  yield  under-estimation  for  suLi-kiloton  weapons  is  more 
severe  when  using  the  expression  from  O'lSA  2506  (Reference  2).  Even 
though  the  author  of  the  report  did  not  provide  an  explicit  caveat 
against  sub-kiloton  application,  we  believe  one  is  implied.  Never- 
theless, the  empirical  equation  of  that  report  is  more  widely  available 
than  the  graphical  results  in  appropriate  classified  reports  and  thus 
frequently  used.  For  range  to  30  psi  peak  overpressure,  the  calculated 
range  is  15  percent  too  large  so  the  yield  would  be  a factor  of  1.5 
too  small.  Modifying  the  analytical  expression  to  use  only  1.65  versus 
a 2.0  multiplier  for  surface  reflection  enhancement,  the  calculated 
range  is  7.5  percent  too  large  for  a yield  factor  of  approximately 
1.25  too  small. 

Another  conclusion  is  that,  for  scaled  range  to  10  psi  peak 
loverpressure,  there  are  statistically  significant  differences  among 
the  BRL,  MOL  and  SC/SRI  instrumentation  types.  For  the  two  higher 
pressures,  instrumentation  type  and  the  interaction  between  yield  and 
instrumentation  type  were  not  significant  contributors  to  the  total 
variance  of  the  data.  (Conclusion  4) 
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The  bias  analysis  of  near-ideal  waveform  (type  V)  versus 
precursor-associated  waveform  types  for  30  psi  peak  overpressure  from 
surface  bursts  did  not  produce  a statistically  significant  result. 
Recognizing  the  subjective  nature  of  the  waveform  parameter,  this 
conclusion  is  also  considered  more  subjective  than  objective  despite 
its  interesting  implications. 
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RANDOM  ERROR  METHODOLOGY 


This  section  describes  the  methodology  developed  to  assess  a 
rationale  for  underground  environment  testing  of  tactical  nuclear  wea- 
pons. based  on  the  assumption  that  variations  observed  during  testing 
are  statistically  the  same  as  variations  that  would  be  experienced  during 
deployment  of  the  weapons. 

The  measure  developed  In  Section  3.1  and  applied  In  later  sections 
to  assess  the  requirements  for  underground  testing  of  tactical  nuclear 
weapons.  Is  the  ratio  of  the  yield  required  to  achieve  a specified  pro- 
bability of  kill  against  a particular  target  assuming  an  Improved  pre- 
dictive capability,  to  the  yield  required  to  achieve  the  same  probability 
of  kill  with  the  present  predictive  capability.  This  measure  Is  developed 
for  both  air  blast  and  nuclear  radiation  effects  on  a point  target  exist- 
ing alone  or  In  a cluster  constituting  a distributed  target.  The  approach 
taken  throughout  attempts  to  obtain  closed-form  expressions  for  the 
measure  when  that  Is  possible;  this  necessitated  making  assumptions  about 
the  distributions  of  weapon  produced  environments  and  target  vulnerability. 
These  assumptions  were  made  on  the  basis  of  reasonableness:  work  has 
been  completed  to  ascertain  the  validity  of  the  assump'’.1on  concerning 
the  distribution  of  the  weapon  produced  blast  environment,  and  Is  re- 
ported In  Section  4.2.  The  development  can,  of  course,  be  generalized 
to  Include  any  assumptions  about  these  distributions,  but  numerical 
solutions  rather  than  closed-form  solutions  might  be  required  or  con- 
venient. 

Section  3.2  describes  the  relationship  between  the  yield  ratio 
measure  and  the  tactical  utility  associated  with  Increased  knowledge 
of  the  nuclear  produced  environment. 


3.1  METHODOLOGY  DEVELOPMENT 


This  Section  builds  the  theoretical  framework  for  assessing 
a rationale  for  underground  testing  of  tactical  nuclear  weapons  based 
on  the  assumption  that  the  weapon  produced  environment  variabilities 
observed  during  testing  are  random  and  will  be  statistically  the  same 
as  the  variabilities  observed  during  weapon  employment.  In  particular, 
expressions  are  derived  for  yield  changes  to  ensure  a given  damage 
level  that  would  result  from  changes  In  the  predictive  capability  for 
weapon  produced  environments. 

To  facilitate  the  development  of  the  expressions  for  yield  ratio, 
a damage  function  Is  derived  which  Is  based  on  the  distributions  for 
the  weapon  created  environment  and  the  target  vulnerability.  This 
function  Is  written  for  convenience  1r,  terms  of  a median  weapon  radius, 
defined  as  the  range  at  which  the  median  weapon  produced  environment 
equals  the  median  target  vulnerability.  The  median  weapon  radius  thus 
becomes  a parameter  of  the  expressions  to  be  developed  for  the  yield 
ratio.  The  damage  function  Is  generalized  to  Include  the  effects  of 
the  weapon  CEP,  and  the  median  weapon  radius  Is  expressed  as  a 
function  of  all  sources  of  uncertainty  (I.e.,  the  uncertainty  In  the 
weapon  produced  environment,  the  uncertainty  In  the  target  vulnera- 
bility, and  the  weapon  CEP)  the  offset  distance,  and  the  kill  proba- 
bility. The  requisite  ratio  of  yields  Is  then  expressed  as  the  ratio 
of  the  above  functions  with  a specified  kill  probability  (Pj^  ■ 0.9), 
where  the  numerator  contains  the  required  uncertainty  In  the  weapon 
produced  environment  and  the  denominator  contains  the  existing  un- 
certainty In  the  weapon  produced  environment.  This  Is  the  measure  used 
to  asses*,  the  effects  of  hypothesized  changes  In  the  uncertainty  of 
weapon  produced  environments  on  the  yield  estimated  to  achieve  a given 
kill  probability.  Each  of  the  following  subsections  treats  some 
pertinent  aspect  of  the  above  development. 
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3.1.1  Derivation  of  the  Damage  Function 


The  damage  function  relates  the  probability  of  damage  at  a given 
range  to  the  distribution  of  the  weapon  produced  environment,  and  the 
distribution  of  target  vulnerability.  For  the  purposes  of  this  develop- 
ment, lognormal  distributions  were  assumed  for  both  the  weapon  produced 
environment  and  the  target  vulnerability.  In  the  case  of  target 
vulnerability,  Reference  1 Indicates  that  for  blast  damage  the  response 
of  the  target  Is  lognormal.  These  data  are  based  on  both  weapons  tests 
and  response  tests  In  controlled  environments.  The  blast  response 
data  have  been  collated  and  reviewed  by  the  Defense  Intelligence  Agency 
and  are  the  basis  of  the  so-called  "VN"  (Vulnerability  Number)  system 
(Ref.  1).  The  radiation  produced  effects  of  Interest  are  prompt  effects; 
the  Information  used  In  this  analysis  comes  from  the  primate  studies 
reported  by  the  Armed  Forces  Radioblology  Research  Institute  (AFRRI) 
(Reference  3). 

In  the  case  of  the  weapon  produced  environments,  no  existing 
source  of  Information  Is  known  to  the  authors  concerning  analysis  to 
Indicate  the  exact  distribution  of  environments.  However,  a lognormal 
assumption  Is  not  unreasonable.  This  hypothesis  Is  addressed  In 
Section  4.2.  Although  the  specific  closed-form  derivations  presented 
herein  depend  upon  assumptions  about  distributions,  the  general 
approach  Is  not,  and  could  be  worked  out  numerically  with  any  choice 
of  distribution. 

To  derive  an  expression  for  the  damage  func  Ion,  let  I repre- 
sent the  random  variable  for  the  nuclear  created  environment  and  V 
represent  the  random  variable  for  the  target  vulnerability,  expressed 
In  the  same  units  as  I.  The  random  variable  I Is  a function  of  ground 
range  R from  the  point  of  detonation,  and  both  the  measure  of  centrality 
of  I (e.g. , mean  or  median)  and  the  variation  of  I (e.g.,  variance) 
could  change  with  R.  Target  damage  occurs  when  I Is  greater  than  V 
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(i.e.,  when  the  weapon  produced  environnvent  Is  greater  than  the  target 
vulnerability).  Defining  a new  random  variable  Z: 

Z - I/V  . (3.1-1) 

The  probability  of  target  damage  at  some  arbitrary  range 


R Is: 


where: 


P(R)  “ f f (Zip)  dZ 


(3.1-2) 


P(R) 


damage  probability 


f (Z|p)  ■ the  probability  density  function  of  Z at  the 
range  R. 

The  Integration  of  f(Zp)is  over  the  range  where  I is  greater  than  V, 
yielding  the  probability  of  target  damage.  If  1 and  V are  lognormally 
distributed  (with  medians  1^  and  V^,  and  variances  s^  and  s^),  then  Z 
Is  lognormally  distributed  with  median  ■ I^/V^  and  variance  s^  « 

*1  * ®v*  terms  of  normal  deviates  the  probability  of  damage 

at  range  R Is  the  result  of  Integrating  Equation  3.1-2. 


rind/z.)-] 

■ f[-^] 


[ 


''q)  1 

s,  J 


(3.1-3) 


A convenient  parameter,  useful  In  the  derivation  of  the  expression 
for  the  yield  ratio  measures.  Is  the  median  weapon  radius,  R^,  defined  as 
that  range  where: 


*0  <«o)  • ''o 


(3.1-4) 
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i.e.,  whe.'e  the  median  weapon  produced  environment  equals  the  median  tar- 
get vulnerability. 

The  next  two  subsections  derive  explicit  expressions  for  the 
damage  function  for  blast  and  radiation  targets. 

3.1.2  Blast  Damage  Function 

For  the  case  of  blast  environments  (peak  overpressure,  overpressure 
impulse,  peak  dynamic  pressure),  the  median  intensity  of  the  weapon  pro- 
duced environment  can  be  assumed  to  follow  an  approximate  power  law, 

Iq  (R)  = Iq  ^*^0^  (R/Ro)‘'‘  • 

for  ranges  close  to  R^.  Thus,  it  is  assumed  that  Iq(R)  can  be 
adequately  represented  by  piecewise  power  laws,  where  the  exponent 
k could  be  different  for  each  section.  Using  Equation  3.1-4,  Equation 
3.1-5  becomes: 

Iq  (R)  = (R/Rq)’*'  =•  Vg  (Rg/R)’'.  (3.1-6) 


The  damage  function  for  blast  targets,  becomes,  on  substituting 
Equation  3.1-6  into  Equation  3.1-3: 


P(R) 


k In  (V^(Rjj/R)/V„) 


(3.1-7) 


where:  E = s^/k,  the  standard  deviation  of  a lognormal  distribution 

R^=  the  median  of  a lognormal  distribution. 

Thus  the  damage  function  in  range  for  blast  targets  is  a lognormal 
distribution  with  variance  E and  median  R„. 
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3.1.3  Radiation  Damage  Function 

For  the  case  of  the  nuclear  radiation  environment,  the  dose 
at  range.  Including  exponential  attenuation  In  air,  can  be  written; 


^0  * ^0  ' 2 In  (R/R^)] 


(3.1-8) 


where  a Is  a range-dependent  Inverse  absorption  length.  Equation  3.1-8 
Is  valid  for  Intervals  around  the  median  range,  R^^,  for  a given  value 
of  a.  By  Equation  3.1-4,  this  equation  becomes: 

Iq  (R)  • V^exp  [a  (R^-R)  - 2 In  (R/R^)]  (3.1-9) 

Using  Equation  3.1-9,  the  damage  function.  Equation  3.1-3,  becomes; 

rin(V  exp[a(R  -R)  - 2 In  (R/R^)]/Vj1 


R„  - R - (2/a)  In  (R/R„)- 


R„  - R - (2/a)  In  (R/R„) 


(3.1-10) 


where:  I » s^/a;  the  standard  deviation  of  a lognormal  distribution 
R|j  ■ the  median  of  a lognormal  deviation 

A good  approximation  to  Equation  3.1-8  for  the  dose  at  ranges  somewhat 
larger  than  the  median  range  (e.g.,  distances  on  the  order  of  1 Km 
for  many  tactical  weapons)  Is  obtained  by  neglecting  the  R-squared  decrease: 

VR)  • loexp(-aR)  ^3  ^ 

This  regime  Is  of  significant  Interest  for  the  tactical  nuclear  case. 

The  radiation  damage  function  assuming  Equation  3.1-11  Is  derived 
analogously  to  the  above,  and  the  result  Is: 

P(R)  - F [-°^|  (3.1-12) 


1 


I 

t 

i 

where  Z is  as  previously  defined  for  Equation  3.1-10;  Equation  3.1-12 
cannot  be  used  near  the  weapon  detonation  point  since  P(R*o)  is  less  than 
one  for  all  R^. 

3.1.4  Derivation  of  Damage  Probability 

The  damage  function  for  point  targets  can  be  generalized  to 
a probability  of  damage  if  offset  targeting  and  the  effects  cf  weapon 
CEP  are  included.  The  case  presented  here  assumes  a circular  bivariate 
normal  distribution  for  weapon  targeting  error  (containing  both 

Urget  anuisition  uncertainties  and  weapon  impact  uncertainties  I 

due  to  guidance,  survey,  etc.).  i! 

Tha  damage  probability  for  a given  CEP,  where  the  CEP  defines  a 
standard  deviation  for  the  bivariate  normal  probability  density  function  | 

(pdf),  where  the  'impoint  is  offset  a distance  d from  the  target  is: 

i 

« -JLj-  Jexp[-{R-d)^/2a^]  P(R)dA  (3.1-13)  | 

2Tfo 

where:  ! 

W * range;  burst  to  target 

P(R)»  damage  function 

- variance  of  the  targeting  pdf;  by  definition  of  the 
CEP  this  variance  is: 

o « CEP//7TnT7r^  0.8493  CEP 

dA  Implies  the  integration  is  ovei  area  (centered  at  the  point 

of  weapon  impact).  i 

Appendix  A presents  the  rather  complex  integration  of  this 
function.  The  results  indicate  that,  for  offset  distance  d > Tr(CEP) 

. p(d)  (3.1-14) 

is  an  adequate  approximation  for  the  damage  probability.  Thus,  for 
offset  distances  greater  than  approximately  three  times  the  CEP,  and 
for  P,j  between  .05  and  .95,  evaluation  of  the  damage  function  at  the 
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offset  distance  is  a valid  approximation  to  Pj^.  For  offset  distances 
<^n  CEP.  Figure  3-1  shows  scaled  median  radius  versus  j (or  versus  damage 
sigma,  Oj)  for  a given  offset,  for  Pj^  =0.9.  In  functional  notation, 
then,  the  scaled  median  radius  can  be  expressed  as: 

yCEP  = g (E,  d/CEP,  P^)  (3.1-15) 

For  Pjj^  ■ .9,  this  function  can  be  written  as: 

g (E,  d/CEP,  .9)  » g (o,  d/CEP,  .9)exp(E^)  (l-G(r))  (3.1-16) 

where  G(l)  is  small  compared  to  one.  Thus,  an  approximation  to  the 
curves  of  Figure  3-1  is: 

Rg  (E)  a Rg(o)  exp(E^)  (3.1-16a) 

or  RJE)*  Rq(o)/(1-o/).  (3.1-16b) 

where  is  the  so-called  damage  sigma  of  the  VN  system  (Reference  1) 
and  is  related  to  I by  the  expression; 


The  next  subsections  utilize  the  above  derived  results  to 
obtain  measures  in  the  form  of  yield  ratios  that  would  indicate  a 
requirement  for  underground  testing  of  tactical  nuclear  weapons  for 
those  cases  where  large  yield  adjustments  could  be  possible. 

3.1.5  Yield  Ratio  Measure 

The  measure  to  be  utilized  to  determine  a rationale  for  under- 
ground testing  of  tactical  nuclear  weapons  is  the  potential  percent 
change  in  yield  by  which  one  would  predict  a given  probability  of 
target  damage  or  a given  level  of  collateral  damage,  if  the  weapon 
produced  environment  were  better  known.  An  alternative  statement  of 
the  measure  is  that  it  is  a measure  of  the  compensation  in  yield  re- 
quired because  the  weapon  produced  environment  is  not  as  well  known 
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Scaled  RMS  Sigma,  Z 
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Figure  3-1.  Scaled  Median  Radius  for  Given  Offset,  , 

P.  “ 0.9,  Blast  Targets  ; 
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as  it  could  be.  An  analytical  formulation  of  this  measure,  applicable 
to  point  targets,  is  the  ratio  of  the  yield,  VI,  required  to  achieve  a 
given  probability  of  damage  (90%)  assuming  the  weapon  produced  environ- 
ment were  perfectly  known,  to  the  yield,  required  to  achieve  the  same 

kill  probability  with  the  present  degree  of  accuracy  in  the  weapon 
produced  environment. 

To  achieve  a significant*  reduction  in  the  yield  predicted 
to  achieve  a 90  percent  kill  probability,  the  ratio  W/W^pj.  must  be  in 
the  range  0.75  to  0.5,  corresponding  to  a 25  percent  to  50  percent 
reduction  in  yield.  The  maximum  improvement  possible  through  UGT  is,  of 
course,  when  the  variance  in  the  weapon  produced  environment  is  set  to 
zero.  To  show  how  the  variance  in  the  weapon  produced  environment  would 
affect  the  yield  selection  to  achieve  a given  probability  of  damage, 
reference  to  Equation  3.1-4  indicates  that  the  median  weapon  radius  is  a 
function  of  yield,  i.e., 

% = f (W)  (3.1-17) 


For  many  weapon  produced  environments  (i.e.,  blast)  R.  is  very  nearly 

1/3  ° 

proportional  to  W ' . Also,  Equation  3.1-15  shows  that  is  a func- 
tion of  I,  the  damage  function  standard  deviation,  where: 

E .=  s^/k  for  blast 

E.  = s^/a  for  radiation 


Thus,  for  those  weapon  produced  environments  where  R is  proportional 

1/3  ® 

to  W , the  ratio  of  yields  becomes: 


EPC 


= (— 
Vd 


-) 


0,  EPC 


rqiz,  d/CEP,  .9) 
'-g(5:^p^,  d/CEP, 


■] 


3 


(3.1-18) 


*Subjective  definition. 
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for  a P|^  of  0.9.  The  value  of  that  will  cause  Equation  3.1-18  to  have 
a value  0.75  o**  smaller  indicates  the  require'  predictive  capability 
potentially  achievable  through  environmental  testing.  Of  course,  there 
may  be  n£  positive  (or  zero)  value  of  s^  tha'  il  result  in  a solu- 
tion in  this  range.  This  implies  one  of  two  iitions  for  UGT: 

t No  significant  gains  in  predictive  capability 
can  be  achieved  through  UGT  if  both  s^  and  ^ 
are  small,  or; 

t Simultaneous  reduction  in  both  s.  and  s is 
indicated,  if  both  and  s^  are  large. 

The  definition  of  "large"  and  "small"  is,  of  course,  subject 
to  interpretation.  Thus,  the  measure  would  clearly  indicate  a 
rationale  for  UGT  if  the  required  significant  yield  change  is  computed, 
but  would  not  necessarily  indicate  no  rationale  if  the  maximum  yield 
change  is  not  significant;  this  later  case  would  require  further 
analysis. 

The  remaining  subsection  derives  explicit  expression  for 
the  ratio  of  yields  for  blast  environment  on  point  targets. 


3.1.6  Yield  Ratio  for  Point  Blast  Targets 

For  point  targets.  Equation  3.1-16  can  be  used  to  express 
the  weapon  radius  as  a function  of Z : 

Rq  (Z)  = Rg  (0)  exp(Z^)  . ' 

The  ratio  of  for  the  requi^-ed  predictive  capability  to  for  the 
existing  predictive  capability  is: 


0,  epc 


= i%  (0)  exp  / [Rjj(O)  exp{Z^p^)] 

= exp[E2  - z^pc] 

5?  + S?  + 

£ exp  [-Lv_  . (_L^_eEc v 


2 2 

s exp  [1J — -iff-] 


)] 


(3.1-19) 


so  that  the  vulnerability  variance, Sy»  cancels  out  of  the  expression. 
Thus,  this  measure  does  not  contain  the  target  vulnerability  when  P|^ 
is  fixed  and  defined  in  terms  of  the  median  weapon  radius. 

A general  expression  for  the  median  weapon  radius  for  blast 
targets  that  encompasses  peak  overpressure  and  dynamic  pressure  and 
overpressure  impulse  is: 


o'  0' 

where:  b * 0 for  peak  overpressure  or  dynamic  pressure, 
b = 1/k  for  impulse. 


(3.1-20) 


Inverting  Equation  3.1-20  to  solve  for  W: 

W 


.A„. 


Thus,  the  ratio  of  yields  becomes: 


(3.1-21) 


epc 


0,  epc 


3 

I+F 


(3.1-22) 
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From  Equation  3.1-19,  the  expression  for  the  ratio  of  median 
weapon  radii  is  inserted  into  Equation  3.1-22: 


n—  = exp 
epc 


[_3_  (ij_ieEc\ 

[l+b  \ ^2  /J  . 


(3.1-23) 


It  can  be  shown  that  the  error  factor  in  intensity,  f,  is 
related  to  the  error  factor  in  range,  fj^,  as: 

f » f!!  . 


(3.1-24) 


Also,  the  error  factor  in  intensity  is  defined: 
f - el-fiS  ^i  . 


(3.1-25) 


Equating  Equations  3.1-24  and  3.1-25  results  in  the  expression  for 


®1  “ TTSS  • 

Substituting  Equation  3.1-26  into  Equation  3.1-23  results  in: 

^R.  epc^]  . 


(3.1-26) 


(3.1-27) 


Equation  3.1-27  is  not  a function  of  target  vulnerability  or 
uncertainty  in  vulnerability.  This  expression  clearly  shows  the 
relationship,  for  all  blast  phenomena  on  point  targets,  between  tl^e 
yield  adjustment  that  one  could  make  to  accomplish  an  objective,  and  the 
Improved  accuracy  in  the  weapon  produced  environment.  The  largest  pos- 
sible adjustment  in  yield  would  be  the  case  where  there  is  no  error  in 
the  estimate  of  the  weapon  produced  environment  (fp,  = 1).  The  value 


of  fg  that  must  ei^fst  so  that  a potentially  significant  reduc- 
tion In^yleld  could  potentially  result  from  environmental  testing  Is 
defined  as  the  critical  value,  rearranging 

Equation  3.1-27: 

^R.  CR  ” 

obtained  by  setting  f^^l  and  Inverting  Equation  3.1-27.  The  significant 
change  In  yield  is  specified  In  Equation  3.1-28  by  setting  W^pc/VJ  equal 
to.  say.  1.25.  indicating  a potential  25  percent  reduction,  and  evalua- 
ting the  expression.  The  result  Is  then  compared  to  fjj,  gp^- 
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3.2 


MEASURE  OF  TACTICAL  UTILITY 


The  tactical  utility  of  a nuclear  weapon  would  be  related  to 
the  area  that  could  be  covered  by  a specified  probability  of  kill,  Pj^, 
given  the  collateral  constraints.  For  any  weapon  there  may  be  a non- 
targetable  excluded  area  defined  by  the  difference  in  the  ranges  between 
the  target-kill  weapon  produced  environment  and  the  col  lateral -damage 
weapon  produced  environment.  If  the  target  is  to  be  killed  with  a 
given  high  probability  (P|^  = 0.9,  for  instance)  and  the  collateral 
weapon  produced  environment  is  to  be  exceeded  only  with  a given  small 
probability  (P„  = 0.1,  for  instance),  then  the  excluded  area  will  be 
a function  of  the  uncertainty  in  these  weapon  produced  environments  — 
knowing  the  weapon  produced  environment  with  a greater  degree  of 
certainty  would  allow  (on  the  average)  a smaller  yield  to  be  specified 
to  achieve  target  kill,  which  would  in  turn  allow  deployment  closer 
to  a collateral  asset,  thus  decreasing  the  excluded  area.  Figure  3-2 
illustrates  the  decrease  in  excluded  area  that  would  result  from 
improved  weapon  produced  environment  predictive  capability. 

To  illustrate  the  relationship  between  this  decrease  in  excluded 
area  and  the  decrease  in  yield  afforded  by  an  improved  predictive 
capability,  consider  the  ratio  of  excluded  areas  for  the  existing  and 
improved  predictive  capabilities: 


where: 

*EPC’  ^IPC  “ 
^EPC’  ^IPC  “ 

•^EPC’  '^IPC  “ 


^EPC 

'iPC 


^EPC  ^EPC 
‘'IPC  '^IPC 


(3.2-1) 


The  areas  for  the  existing  and  improved  predictive 
capabilities. 

The  ranges  to  the  collateral  damage  weapon  produced 
environment  for  the  existing  and  improved  pre- 
dictive capabilities. 

The  ranges  to  the  target  kill  weapon  produced 
environment  for  the  existing  and  improved  predictive 
capabilities. 
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The  blast  environment  at  a given  range  close  to  the  median  weapon  radius 
was  given  by  Equation  3.1-6  as: 


i(R)  = . 

(3.1-6) 

Let  I be  some  specified  weapon  produced  blast  environment. 
Equation  3.1-6  for  the  range  to  this  specified  environment 

Solving 

results  In: 

R,  ■ Ro(v„/i)‘''‘  . 

(3.2-2) 

Thus  the  ranges  to  specified  target  and  collateral  environments  are, 
respectively: 

R ■ Ro(V‘t>‘''‘ 

(3.2-3) 

c - RqCVc/Ic)^^''  . 

(3.2-4) 

Equation  3.1- 16a  gives  a conceptual  form  for  In  terms  of  the  weapon 
CEP,  the  specified  Pj^,  and  the  weapon  and  target  vulnerability 
uncertainties  as: 

Rq  “ g{o,  d/CEP,  P,^)e^  . (3.1-16a) 

Actually  this  expression  Is  valid  only  for  Pj^'s  near  0.9;  a more 
general  expression,  encompassing  both  target  kill  and  collateral  damage 
P|^'s,  would  be: 

Rq  - g(o,  d/CEP,  P,^)f{Z)  (3.2-5) 

where  f(2)  would  be  different  for  target  and  collateral  P|^'s.  On 
substituting  Equation  3.2-5  Into  Equations  3.2-3  and  3.2-4,  they 
become : 
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\ , 
■ \ 


(3.2-6) 


R - 9(0.  d/CEP. 

C - g(o,  d/CEP.  * 


(3.2-7) 


Equations  3.2-6  and  3.2-7  can  now  be  used  In  Equation  3.2-1  to  express 
the  ratio  of  excluded  areas  In  terms  of  the  target  and  collateral  Pj^'s 
and  the  existing  and  Improved  predictive  capabilities.  Equation  3.2-1 
becomes: 


'VA2/k 


‘EPC 

iPC 


g2(o.  d/CEP.  Pk)f^J:EPc(T;) 

» 1'^ 

9^(0.  d/CEP.  Pc)fc(^iPcHT 


9^(0.  d/CEP.  Pk)f^(%cHri 

r 

g2(o.  d/CEP, 


(3.2-8) 


cut  the  range  to  a particular  nominal  effect  Is  related  to  the  yield 
as: 


R " h W 


,1/3 


(3.2-9) 


Where  k Is  the  range  to  the  effect  for  blast.  Substituting 
Equation  3.2-9  Into  Equation  3.2-8,  where  the  yields  Wgpj;  and  Wjpc 
have  the  relationship  described  In  Section  3.1,  results  In: 


'V„\2/k 


-EPC 

'iPC 


k2  w 2/3 
"c  *‘lPC 


,2/3 


rxm 


(3.2-10) 
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Thus,  the  ratio  of  the  excluded  areas  for  the  existing  and  improved 
predictive  capabilities  is  equal  to  the  two-thirds  power  of  the 
yield  ratio. 
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Section  4 

RANDOM  ERROR  RESULTS 


This  section  describes  the  results  of  the  assessment  of  a 
rationale  for  tactical  UGT  environment  tests  based  on  the  assumption 
that  the  nuclear  produced  blast  environments  observed  in  the  historical 
test  data  are  random,  and  are  statistically  the  same  as  the  environ- 
ments to  be  expected  on  deployment  of  tactical  weapons.  Section  4.1 
presents  the  blast  data,  and  the  results  for  the  rationale  assessment 
based  on  analysis  of  peak  oven.ressure  data,  overpressure  impulse  data, 
and  peak  dynamic  pressure  data.  Section  4.2  presents  the  results  of 
the  investigation  of  the  statistical  distribution  appropriate  for 
these  data. 

4.1  RANDOM  ERROR  ANALYSIS  OF  BLAST  DATA 

This  section  presents  the  raw  blast  data  and  shows  the  results 
of  applying  the  measure  developed  in  Section  3 to  assess  a rationale 
for  underground  nuclear  effects  tests  This  methodology  was  developed 
in  Section  3 in  closed  form  under  the  assumption  that  errors  in  the 
scaled  range  to  a specified  environment  are  proportional  to  range. 
Section  4.2  examines  the  validity  of  that  assumption  in  some  detail 
(the  conclusion  is  that  the  assumption  of  lognormality  is  a reasonable 
and  adequate  one  for  most  nuclear  blast  data). 

Figures  4.1-1  through  4.1-31  show  the  peak  overpressure,  over- 
pressure Impulse,  and  peak  dynamic  pressure  data  plotted  versus  scaled 
range  for  all  ten  ranges  of  scaled  HOB.  These  plots  show  the  scatter 
in  the  "raw"  data,  as  well  as  the  approximate  linear  relationship  when 
plotted  on  a log-log  scaled  coordinate  system,  illustrating  the  approx- 
imate power-law  relationship  between  these  three  blast  effects  and 
scaled  range. 

The  measure  previously  developed  applies  to  the  scatter  in  the 
range  to  a specified  environment,  and  is  based  on  the  range  error 
factor  at  the  specified  environment.  To  compute  the  range  error  factors 
at  specified  environments,  it  is  necessary  to  group  the  blast  data  into 
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Figure  4.1-1.  Peak  Overpressures  for  0 to  3 ft/kt^'^^  HOB 
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Figure  4.1-4.  Peak  Overpressures  for  113  to  157  ft/kt^'^^  HOB 
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Figure  4.1-5.  Peak  Overpressures  for  182  to  205  ft/kt^^^  HOB 
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Figure  4.1-6.  Peak  Overpressures  for  212  to  252  ft/kt^^^  HOB 
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Figure  4.1-7,  Peak  Overpressures  for  323  to  375  ft/kt^^^  HOB 
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Figure  4.1-8.  Peak  Overpressures  for  478  to  500  ft/kt^^^  HOB 
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Figure  4.1-11.  Overpressure  Impulses  for  0 to  3 ft/kt^^"^  HOB 


BEST  AVAIIABIE  COPY 


Overpressure  Impulse  (psl-msec/kt^^^) 


Overpressure  Impulse  (ps1-msec/kt^^^ 


Ground  Range  (ft/kt^^^) 

Figure  4.1-13.  Overpressure  Impulses  for  55  to  83  ft/kt^*^^  HOB 
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Figure  4.1-14.  Overpressure  Impulses  for  113  to  157  ft/kt^^^  HOB 
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Figure  4.1-15.  Overpressure  Impulses  for  182  to  205  ft/kt^^^  HOB 
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Figure  4.1-i7.  Overpressure  Impulses  for  323  to  375  ft/ HOB 
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Figure  4.1-19.  Overpressure  Impulses  for  751  to  831  ft/ktM^  HOB 
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Figure  4.1-20.  Overpressure  Impulses  for  1003  to  1249  ft/kt^^^  HOB 
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Figure  4.1-21.  Peak  Dynamic  Pressure  for  137  and  156  ft/kt^'^^  HOB 
(Heavy  Dust) 
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Figure  4.1-22.  Peak  Dynamic  Pressure  for  188  and  204  ft/kt^^^  HOB 
(Heavy  Dust) 
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Figure  4.1-24.  Peak  Dynamic  Pressure  for  0 to  2.7  ft/kt^^^  HOB 
(Light  Oust) 

4-25 


best  AVAIIABIE  COPV 


iiikiAifiMri’  I 


Peak  Dynamic  Pressure  (psi) 


40& 


160i 

I 

looi 


40  V” 

loi 

I 

« 

4 

t 

• 

I 

I 

4l  — 

i.el 

i.ol 

I 

4 

• 

0.4L_ 

200 


a 


Figure  4.1-25. 
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Figure  4.1-26.  Peak  Dynamic  Pressure  for  81  and  113  HOB 
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Figure  4.1-28.  Peak  Dynamic  Pressure  for  180  to  196  ft/kt"'"'  H08 
(Light  Dust) 
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Figure  4,1-31.  Peak  Dynamic  Pressu."  for  786  ft/kt^^^  HOB 
(Light  Dust) 
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range  bins.  Both  coarse  and  fine  range  bins  were  chosen  for  this 
analysis  so  that  any  effects  that  could  result  from  a particular 
grouping  of  these  data  would  be  illustrated  (none  were  found;  see 
Section  4.2).  Bin  sizes  corresponding  to  12  bins  per  decade  and  3 
bins  per  decade  were  chosen  to  group  these  data.  To  eliminate  trends 
within  bins,  due  to  the  power  law  relationship  between  range  and  the 
nuclearly  produced  environment  within  the  bin,  the  data  within  each 
bin  were  adjusted  by  the  approximate  power  law  relationship. 

Tables  4.1-1  through  4.1-3  present  the  range  error  factors 
computed  for  each  bin  for  both  peak  pressures  and  overpressure  impulse 
and  for  the  scaled  HOB  ranges  shown  on  the  tables.  Tables  4.1-1  and 
4.1-2  show  the  range  error  factors  computed  by  grouping  these  data 
into  12  bins  per  decade  for  peak  overpressure  and  overpressure  impulse, 
respectively,  and  Table  4.1-3  presents  error  factors  for  both  of  the 
above  nuclear  produced  environments  and  peak  dynamic  pressure  computed 
by  grouping  these  data  into  3 bins  per  decade.  In  addition,  the  tables 
show  the  range  error  factor  averaged  across  SHOB's  for  each  range  bin 
(row  averages)  ana  across  range  bins  for  each  SIlOB  bin  (column  averages). 
These  error  factors  represent  the  variation  in  these  data  upon  which 
the  yield  ratio  measure  developed  in  Section  3 was  applied. 

To  compare  the  results  for  the  different  size  bins,  range 
error  factors  for  each  bin  size  were  plotted  versus  range,  for  each 
of  the  scaled  HOB's  in  the  tables.  Figures  4.1-32  through  4.1-41 
show,  for  each  scaled  HOB  bin,  the  superimposed  plots  of  range  error 
factor  versus  range  for  the  peak  overpressure  data  for  both  bin  sizes. 
Figure  4.1-42  through  4.1-51  are  similar  plots  for  overpressure 
impulses,  and  Figures  4.1-52  and  53  show  range  error  factors  versus 
range  for  peak  dynamic  pressure  (combining  these  light  and  heavy  dust 
data  would  be  possible  since  there  appears  to  be  no  difference  in 
range  error  factors).  These  data  are  summarized  in  Figure  4.1-54, 
which  shows  the  average  peak  overpressure  and  overpressure  impulse 
range  error  factors  versus  range  for  both  types  of  binning.  Figures 
4.1-32  through  4.1-54  indicate  the  following: 
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t The  grouping  of  these  data  into  bins  of  various  widths 
does  not  appear  to  appreciably  affect  the  estimates 
of  range  error  factors,  since,  in  practically  all  cases 
the  3-bin  per  decade  and  12-bin  per  decade  curves  are 
coincident. 

• For  peak  overpressure  and  dynamic  pressure,  the 
plots  of  range  error  factor  versus  range  appear  to 

be  constant.  For  overpressure  impulse,  there  appears 
to  be  a slight  decrease  in  range  error  factor  with 
increasing  range. 

e The  variation  in  range  error  factor  from  bin  to  bin, 
and  across  scaled  HOB  is  not  large;  there  appear  to 
be  few  anomalous  cases. 

All  of  the  above  observations  would  suggest  that  there  appear  to  be 
few  "special  cases"  where  the  yield  ratio  measure  would  indicate  a 
substantially  different  conclusion  from  that  based  on  the  measure 
applied  to  the  average  data. 

The  minimum  and  maximum  range  error  factor  values  corresponding 
to  the  following  combinations  of  characteristics  in  Tables  4.1-1 
through  4.1-3  were  chosen  to  compute  the  yield  ratio  measure: 

• The  cell  minimum  and  maxiumum  for  peak  overpressure, 
overpressure  impulse,  and  peak  dynamic  pressure  for 
3 and  12  bins  per  decade. 

t The  minimum  and  maximum  averaged  for  each  column 
(corresponding  to  a particular  SHOB)  for  peak 
overpressure,  overpressure  impulse,  and  peak  dynamic 
pressure,  for  both  types  of  binning. 

• The  minimum  and  maximum  averaged  for  each  row 
(corresponding  to  a particular  range  bin)  for  peak 
overpressure,  overpressure  impulse,  and  peak  dynamic 
pressure,  for  both  types  of  binning. 

Table  4.1-4  presents  both  the  range  error  factor  and  the  yield  ratio 
measure  for  these  characteristic  data.  The  yield  ratio  measures 
larger  than  10  percent  (i.e.,  that  would  suggest  a possible  maximum 
10  percent  or  larger  reduction  in  yield  through  additional  under- 
ground environments  testing)  are  displayed  in  Table  4.1-5  along  with 
a brief  discussion  of  the  possible  tactical  significance  of  effecting 
a yield  reduction. 
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Table  4.1-1  RAflGE  ERROR  FACTORS  PEAK  OVERPRESSURE  (12  BIHS/DECADE) 


o 

o 


Table  4.1-2  RANGE  ERROR  FACTORS  OVERPRESSURE  IIIPULSE  (12  BINS/OECADE) 


RANGE  ERROR  FACTORS  - PEAK  OVERPRESSURE,  OVERPRESSURE 
IMPULSE.  AND  PEAK  DYNAMIC  PRESSURE  (3  BINS/OECAOE) 


3 bins/decade 
2 bins/decade 


■ 3 bins/decade 
12  bins/decade 
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Figure  4.1-50.  Range  Error  Factor  vs.  Scaled  Range 

SHOB:  751-831  ft  (Overoressure  Impulse) 
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SHOB: 


Table  4.1-4  SUWARY  OF  YIELD  RATIO  MEASURE  FOR  81  AST  ENVIRONMENTS 
(by  cell,  SHOB.  and  rangt  bin) 
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bad"  SUGAR  datum  having  large  Impact  on  the  error  factor  was  ignored 


Table  4.1-5  IMPACT  OF  ANALYSIS  ON  TACTICAL  UTILITY 
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DISTRIBUTION  TYPE  FOR  BLAST  DATA 


4.2.1  Introduction 

The  measurements  of  scaled  range  to  a selected  pressure  value 
have  a distribution  which  can  be  statistically  described.  The 
distribution  must  be  characterized  in  order  to  calculate  the  scaled 
range,  to  the  selected  pressure  value,  with  high  and  low  probabilities 
for  target  and  collateral  damage  analyses  respectively.  Since  the 
total  population  is  required  to  specify  the  exact  distribution,  the 
existing  sample  can  only  be  used  to  specify  a best  distribution. 

The  following  analysis  was  performed  to  find  the  dependence 
of  range  uncertainty  on  range.  The  hypothesis  is  that  a linear  or 
higher  order  dependence  completely  negates  the  applicability  of  a 
normal  distribution.  Specifically  a lognormal  distribution  applies 
for  a linear  dependence.  For  higher  order  dependence,  a linear 
dependence  approximation  is  usually  dominant  and  thus  assumed  to 
apply,  but  not  to  the  total  exclusion  of  another  continuous  non- 
normal  distribution  such  as  beta  and  gamma  which  have  non-zero 
skewness.  The  lognormal  is  more  convenient  than  beta  or  gamma 
because  the  data  are  frequently  plotted  on  a logarithmic  scale 
and  the  data  uncertainty  is  then  a symmetric  factor  in  linear  space. 

The  measurements  used  in  this  analysis  are  the  scaled  range 
to  static  (over)  pressure  peak  and  scaled  impulse  values.  A data 
quality  control  was  attempted  by  selecting  only  the  ranges  where  both 
peak  and  impulse  values  were  simultaneously  measured.  The  data  base 
is  that  previously  culled  set  which  is  reported  in  DASA  1200,  Volume  V 
(Reference  4).  In  addition,  the  biased  data  (Reference  D)  of  the 
Franklin  event  was  not  used  The  similar  data  base  for  dynamic 
pressure  is  unpublished  but  has  been  partially  compiled  by 
Mr.  J.  Keefer,  et  al.,  of  the  U.  S.  Army  Ballistics  Research  Labora- 
tory. This  dynamic  set  is  nuch  smaller  than  the  static  set  so  esti- 
mating the  standard  deviation  would  not  be  as  accurate  and  thus  the 
hypothesis  p'^oof  is  more  conn  lex. 


The  approach  for  calculating  the  standard  deviation  of  scaled 
range  as  a function  of  scaled  range  is  the  following.  The  measurements 
are  systematically  grouped  or  binned  in  scaled  range.  The  data 
within  a range  bin  are  normalized  to  the  range  at  the  mid-point  of 
the  bin  using  a power  law  approximation.  Then  these  range-square- 
normalized  peak  overpressures  and  range-normalized  overpressure 
impulses  are  analyzed  for  a mean  and  variance.  This  overpressure 
variance  is  then  propagated  to  a range  variance  using  the  law  of 
covariance  propagation; 


- im 
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where  z is  either  overpressure  peak  or  impulse. 


The  second  term  involving  the  third  central  moment  (skewness 
indicator)  can  be  ignored.  The  few  cases  having  an  adequate  number 
of  values  (i.e.,?  15)  in  a bin  to  estimate  skewness  have  both 
positive  and  negative  contributions.  The  standard  deviation  change 
attributable  to  the  correction  for  skewness  was  less  than  15  percent 
which  can  later  be  inferred  as  insignificant. 


The  final  steps  consist  of  plotting  the  standard  deviations 
of  range  versus  the  mid- points  of  the  range  bins  and  fitting  a linear 
equation  through  the  Vflues.  The  uncertainty  of  the  slope  coefficient 
allows  testing  for  significance  of  the  linear  fit.  The  following 
material  discusses  the  details  of  this  procedure  and  the  results. 


4,2.2  Procedure 


The  lack  of  a large  amount  of  data  for  a single  scaled  HOB 
other  than  Priscilla  requires  SHOB  grouping  or  binning  to  accumulate 
a large  number  of  data  within  each  range  bin.  Since  a simple  SHOB 
normalization  for  all  types  of  blast  data  is  unavailable,  the  SHOB 
bins  are  kept  small.  The  data  base  was  compartmented  into  10  SHOB 
bins  using  the  bar  graph  of  Figures  4.2-1  through  4.2-4  as  guides. 


4-64 


Scaled  HOB 


The  SHOE  range,  number  of  simultaneous  peaks  and  impulses,  and  names 
of  tests  in  the  10  bins  are  as  follows: 


HOB  - ft/kt^^^{ID) 

No.  of 
Data 

Events 

0 - 3.1  (03) 

177 

Mike,  Walnut,  Yankee, 

Union,  Romeo,  Zuni, 

Koa,  Bravo,  Nectar 

Cactus,  Yellowwood,  Aspen, 
Sycamore,  Koon,  Sugar 

5.0  - 10.8  (511) 

87 

Lacrosse,  Fig,  Small  Boy, 

Little  Feller  I & II 

55  - 83  (5583) 

104 

Oog,  Zebra,  Simon,  Inca, 

Easy,  Yoke,  X-Ray 

113-157  (113157) 

62 

Annie,  Turk,  Met,  Humboldt, 
Zucchini,  Tesla,  Apple  2 

182  - 205  (182205) 

143 

King,  Shasta,  Hornet,  Smoky, 
Apple  1,  Priscilla,  Grable 

212  - 252  (212252) 

96 

Galileo,  Moth,  Kepler,  Wilson, 
Owens,  Morgan,  Bee,  Post 

323  - 375  (323375) 

81 

Climax,  Hood,  Yuma,  Dog 

478  - 500  (478500) 

41 

Wasp  Prime,  Ruth,  Hamilton 

751  - 831  (751831) 

65 

Soccorro,  Able,  Encore,  Rushmore 

1003  - 1249  (10031249) 

97 

Baker,  Charlie,  Eddy,  Mora,  Lea 

Specific  events  with  numerous  data  in  the  1 to  100  psi  interval  are 
the  following: 


1.2 

28 

Cactus 

8.2 

36 

Small  Boy 

83 

37 

Easy 

204 

57 

Priscilla 

375 

40 

Dog 
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The  data  measured  near  the  surface  are  for  both  regular  and 

Kach  stem  reflections  of  the  initially  spherical  air  blast  wave.  The 

knee  at  the  transition  between  the  reflection  regions  on  a olot  with 

HOB  and  ground  range  coordinates  is  very  nonlinear  and  not  directly 

coirtjatible  with  a statistical  treatment.  Most  targeting  analyses  use 

the  Mach  stem  overpressure  enhancement  so  only  the  data  below  the  knee 

are  used,  in  this  tactical  analysis,  for  air  burst  events.  The  bottom 

of  the  knee  used  to  isolate  these  data  was  described  by  the  approximate 

1/3 

empirical  equation  for  units  of  ft/kt  : 

0 2000  X IlOB 

“ ■ 200Cr-~il0B 

where  R is  222  for  a 2G0  IlOB  and  R is  500  for  a 400  IlOB. 

The  variance  has  both  a true  and  an  error  component  unless  all 
points  have  the  same  mean.  This  same  mean  is  achieved  by  removing 
the  analytical  dependence  on  range  in  addition  to  HOB.  The  error 
variance  is  very  small  by  using  a power  law  approximation  for  normali- 
zation within  the  range  bin.  Peak  overpressures  were  normalized  by 
inverse  range  squared  and  overpressure  impulses  were  normalized  by 
inverse  range.  This  approximate  power  of  two  for  peak  overpressure 
could  be  refined  by  using  the  Erode  Direct  Fit  (Reference  2)  to  estimate 
the  power.  But  the  Brode  pressure  time  equations  are  not  adequately 
suitable  for  impulse  calculations  where  the  scaled  measurements  were 
made  so  range  normalization  of  impulses  would  not  benefit  by  using  Brode's 
Direct  Fit. 

Another  more  direct  approach  to  minimizing  the  error  variance 
is  to  use  very  narrow  range  bins  so  the  range  normal ization  is  negligible. 
This  also  has  the  effect  of  emphasizing  duplicate  measurements  at  the 
same  point  which  provide  a true  variance.  Both  linear  and  logarithmic 
range  bins  were  used  to  force  different  groupings  of  the  data  as  the 
range  bin  size  was  varied.  The  sizes  for  both  types  of  bins  were 
varied  by  a factor  of  ten  between  largest  and  smallest  with  minor 
differences  in  the  calculated  variances  and  no  apparent  biases  or  trends. 


The  square  roots  of  some  variances  for  group  iSCS  are  plotted  in  Figure 
4.2-5  and  isted  in  Table  4.2-1.  Tnis  result  consistency  is  verifica- 
tion of  the  adequacy  of  the  simple  pov.-er  lav;  range  normalizations  for 
peak  and  iiTDulse  values  within  each  bin. 

The  possible  error  variance  from  nixing  data  from  precursor 
and  non-precursor  measurements  is  also  of  interest.  For  this  immediate 
analysis,  the  interest  is  only  in  the  possible  systematic  bias  of  the 
variance  with  range.  Since  neither  the  respective  means  nor  thair 
difference  are  precisely  defined,  the  issue  was  indirectly  treated  by 
assessing  the  representativeness  of  the  true  variance  from  a given 
event  by  the  variance  front  the  SllOB  group  including  that  event.  Spec- 
ifically, event  Easy  from  group  55T3  and  event  ^’riscilla  from  group 
182205  were  used  for  this  intra-comoar ison.  The  standard  deviations 
of  peak  overpressures  are  plotted  in  Figures  4.1‘-6  and  4,2-7  co  show 
the  representativeness,  ‘iote  chat  some  event  variances  are  larger 
than  the  respective  group  variances  and  no  bias  or  trends  are  notice- 
able. Too  few  values  exist  for  a linear  fit  ev  the  event's  variances 
to  allow  reliable  quantitative  cor.iparisor  ■ ith  a linear  fii  of  the 
groun's  variances. 

fhe  sparsity  of  data  causes  concern  for  the  reliability  of 
variance  estimates.  Five  points  were  chosen  (approximately  one  in  the 
middle  and  two  on  both  sides)  to  help  maximize  the  number  of  variance 
estimates.  Using  the  variance  as  the  accepted  measure  of  data  dis- 
persion, not  as  the  estimator  of  the  standard  deviation  (which  implies 
normality),  reduces  the  concern  for  mixing  variances  from  five  and, 
for  example,  twenty-seven  points.  Four  points  are  expected  to  all  be 
within  one  standard  deviation  and  twenty-seven  points  are  expected  to 
all  be  within  two  standard  deviations  for  points  normally  distributed 
around  an  average.  This  guidance  provides  analytical  justification  for 
dropping  points  which  differ  drastically  from  the  average,  that  is, 
outliers.  Empirical  justification  for  dropping  outliers  was  also 
necessary  in  the  form  of  severe  departure  from  the  trend  line  of  the 
other  points  for  the  same  event. 
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Deviation  of  Peak  Overpressure 


Table  4.2-1.  Standard  Deviations  of  Peak  Overpressures  and  Overpressure  Impul 
for  Linear  and  Logarithmic  Range  Bins 
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Standard  Deviation  of  Peak  Overpressure  (psi) 
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The  representativeness  of  as  few  as  five  points  can  also  be 
discussed  without  reference  to  any  assumption  about  the  type  of  distri- 
bution (Reference  6).  The  result  is  expressed  as  f'.e  continuous  prob- 
ability, with  nur.ber  of  points  in  the  sample  as  parameter,  for  a frac- 
tion of  the  population  being  within  the  sample  extremes.  The  expected 
value  of  the  fraction,  from  many  samples  each  of  n points,  is  (i-l)/ 
(n+1)  but  the  probability  of  this  value  of  the  fraction  slowly  increases 
from  0.5  for  n=3  to  an  asymptotic  value  below  0.6.  Thus  the  0.542 
probability  of  including  67  percent  of  the  population  with  a sample  of 
5 points  is  not  too  different  than  the  0.536  probability  of  including 
93  percent  of  the  population  with  a sample  of  27  points.  See  the 
expected  value  curve  in  Figure  4.2-3.  Also  plotted  on  the  graph  is 
the  more  interesting  dashed  curve  labelled  most  probable.  The  most 
probable  value  of  the  fraction,  from  one  sample  of  n points,  is  (n-2) 
/(n-1),  but  the  probability  of  this  value  of  the  fraction  steadily 
decreases.  Thus  one  is  not  automatically  assured  of  having  a fraction 
larger  than  that  for  the  expected  value.  The  net  result  is  that  5 
points  was  selected  as  a minimum  which  implies  a minimum  of  4 degrees 
of  freedom  for  the  variance  estimate. 

The  bin  centers  were  also  staggered  to  maximize  the  number  of 
variances.  Two  runs  of  an  identical  bin  size  were  made  with  50  percent 
offset  of  each  bin  with  respect  to  the  bins  in  the  other  run.  This 
effort  produced  an  algorithm  which  could  successfully  bin  and  calculate 
variances  from  more  than  90  percent  of  the  data.  The  results  for  group 
5583  are  tabulated  as  follows: 


Bins/Decade 

46 

23 

15 

12 

8 

5 

3 


No.  of  Max  No. /Bin 


8 

11 

10 

7 

6 

4 

3 


10 

16 

20 

23 

30 

38 

52 


No.  of 

11 

12 

9 

1 8 
5 
4 
2 


Max.  No/Bin 

10 

12 

16 

23 

27 

41 

62 
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Probability  of  a for  Sample  Size  N 


a Fraction  of  Pooulation  Within  Samole  Extremes 


Figure  4.2-8.  Probability  of  a Fraction  of  the  Population  Being 
Within  Sample  Extremes 
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Note  that  even  the  4 variances  are  sufficient  to  least  squares  fit  a 
linear  curve  and  test  for  a statistically  significant  slope. 

This  group  of  104  points  has  some  range  redundancy  as  evident 
by  the  large  number  of  bins  with  5 through  10  points  for  46  bins  per 
decade  (that  is,  the  ratio  of  maximum  to  minimum  ranges  for  the  bin  is 
only  1.051).  This  much  data  prevents  the  offset  feature  from  operating 
effectively.  For  event  Easy  (only  37  points)  at  the  maximum  SHOB  of 


this  group,  the  following 

tabulation  results: 

Bins/Decade 

No.  of 

Max.  No. /Bin 

No.  of 

Max.  No. 

46 

1 

8 

1 

8 

23 

2 

10 

2 

8 

15 

3 

14 

3 

8 

12 

2 

10 

4 

14 

8 

2 

14 

4 

14 

5 

3 

14 

3 

16 

3 

2 

20 

2 

28 

Thus  the  final  algorithm  is  applicable  to  many  individual  events 
where  independent  blast  instrumentation  lines  were  used.  An  empirical 
estimate  of  the  maximum  number  of  variances  from  these  two  examples  is 
one-ninth  of  the  number  of  points. 

The  next  step  required  obtaining  an  approximate  analytical 
relationship  for  differentiating  in  order  to  convert  pressure  variance 
to  range  variance.  The  log  p - log  scaled  R data  was  linearly  regressed 
with  attention  on  the  coefficient  of  determination  (squared  regression 
coefficient)  which  measures  the  amount  of  variation  from  the  mean 
explained  by  the  fitted  curve.  (Alternatively,  a "forecast  efficiency" 
is  calculated  from  1 - (1-r^)*^  which  is  plotted  in  Figure  4.2-9;  only 
50  percent  efficiency  results  from  r = 0.866  or  r = 0.75.)  Note  also 
that  the  coefficient  of  determination  is  the  product  of  the  slopes  from 
independently  regressing  y-on-x  on  x-on-y  so  it  implies  the  inaccuracy 
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of  using  the  inverse  of  the  p-on-R  (i.e.,  Erode  equacion)  for  a R-on-p 
•^it  which  is  needed  for  differentiation  and  use  in  the  variance 
oropagation.  Recognizing  that  one  straight  line  across  all  the  data 
is  wrong,  the  linear  regression  was  performed  in  independent  intervals 
using  very  wide  bins;  a quadratic  fit  is  more  appropriate  but  causes  a 
messy  differentiation.  The  results  for  various  data  and  bin  widths 
follow  in  Table  4.2-2. 

Table  4.2-2.  Coefficient  of  Determination  for 
Broad  Range  Bins 


2 Bins/Decade 

3 Bins/Decade 

5 Bins/Decade 

*^med 

N rj 

Ap 

•"ip 

*^med 

N 

2 2 
"ap  "ip 

*^med 

N 

i 

Group  5585 

1 

562 

66  .86 

.51 

316 

14 

.28  .01 

398 

17 

.15 

.01 

1778 

40  .74 

.63 

681 

52 

.79  .42 

631 

33 

.66 

.13 

1469 

40 

.74  .63  1 

1000 

41 

.35 

.30 

1585 

15 

.87 

.61 

Priscilla 

278 

6 .62 

.42 

316 

29 

.96  .59 

316 

22 

.93 

.50 

562 

46  .96 

.64 

681 

23 

.82  .43 

501 

21 

.81 

.10 

1778 

5 .39 

.04 

1469 

5 

.39  .04 

794 

9 

.56 

.31 

1 

1259 

5 

.39 

.04 

j 

Group  212252 

562 

63  .88 

.82 

464 

29 

.94  .81 

316 

6 

.92 

.86 

1778 

28  .92 

.93 

1000 

54 

.47  .67 

501 

18 

.73 

.59 

2153 

8 

.98  .96 

794 

39 

.16 

.24 

1259 

21 

.67 

.58 

1995 

6 

.97 

.89 
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These  results  have  a clear  message.  The  peak  overpressure 
measurements  are  much  better  described  by  a power  law  than  are  the 
overpressure  impulse  measurements.  The  consequence  is  that  the 
following  quadratic  form  should  be  used  for  overpressure  impulse  without 
range  binning  for  all  the  data: 

For  peak  overpressure,  the  self-consistency  of  the  data  is  not  adequate 
to  prefer  these  fits  to  the  Brode  equation  as  far  as  calculating  the 
derivative  or  slope  for  t le  variance  propagation.  Thus  the  two  peak 
overpressures  were  calci  ited  from  the  Brode  equation  for  bin  edges 
and  the  finite  different  approximation  calculated  for  the  slope  at 
the  mid-point.  The  second  derivative  was  likewise  approximated  so  the 
ratio  of  the  skewness  term  to  the  variance  term  could  be  calculated 
and  used  for  propagation  of  peak  overpressure  scatter  to  scaled  range 
variance. 

The  scaled  overpressure  impulse  propagating  in  free  air  is 
inversely  proportional  to  scaled  range  (Reference  7).  and  this  simple 
power  law  approximation  was  adequate  for  range  normalization  within 
both  linear  and  logarithmic  bins  (see  Figure  4.2-10  for  consistency  of 
variance  estimates  from  narrow  range  bins).  Thus  the  inverse  linear 
power  is  assumed  so  the  derivative  is  as  follows; 


which  leads  to  the  propagation  formula; 


a 


RI 


or  Ca^j 


4-81 


- / 

/ 


standard  Deviation  of  Overpressure  Impulse  (psi-msec/kt^/^) 


The  proportionality  constant  C becomes  a scale  factor  and  can  be  ignored 
later. 

The  final  step  is  the  fitting  of  the  range  uncertainty  to  range. 
The  highest  correlation  is  expected  by  fitting  in  logarithmic  space. 

But  the  lognormal  hypothesis  requires  fitting  specifically  in  linear 
space.  That  is,  the  polynomial 

= a + bR 

is  used  and  b is  tested  for  significance  using  its  standard  deviation 
calculated  from  the  remaining  average  residual. 

4.2.3  Results 

Most  of  the  data  base  of  static  pressures  is  presented 
graphically  in  the  20  plots  of  Section  4.1  for  various  scaled  HOB  bins. 
Only  peaks  with  associated  impulses  were  used  so  the  results  could  be 
intercompared  directly.  The  sparce  amount  of  data  beyond  the  edges  of 
the  plots  were  neglected  only  for  plotting  convenience.  The  minimum 
and  maximum  axes  values  cover  the  region  of  t'ctical  interest,  but  the 
scientific  interest  extends  another  factor  of  ten  or  more. 

This  analysis  has  emphasized  data  in  Reference  4 which  was 

further  culled  so  systematic  processing  could  be  used.  The  culled  data 

were  primarily  impulses.  The  following  tabulation  presents  one  of  the 

author's  (JEC)  subjective  reasons  for  culling  the  value.  A small 

1/1  ” 
value  of  145  at  2064  ft/kt  ' for  Sugar  was  not  culled  so  its  effect 

can  be  followed  through  the  computerized  analysis. 


Series 

Event 

R(ft/kt^^^) 

Subjective  Reason 

Hardtack  1 

Koa 

296 

Both  Ip  too  large 

Redwing 

Zuni 

454 

Ip  too  small 

Castle 

Bravo 

521 

Ip  too  small 

Hardtack  1 

Fig 

253,  326 

Ip  too  large 

Sunbeam 

Small  Boy 

313 

Ip  too  large 

Greenhouse 

Easy 

383,  466 

1^  too  large 

Tumbler 

Dog 

794 

Ip  too  large 

Hardtack  II 

Hamilton 

1839 

Ip  too  small 

The  standard  deviations  of  range  were  normalized  to  the  smallest 
one  in  the  respective  SHOB  group  and  this  ratio  is  plotted  in  the  next 
set  of  nine  graphs.  Figures  4.2-11  through  4.2-19.  (Too  fev  Jata  exist 
in  this  analysis  of  SHOB  group  478500  for  the  8 or  12  bins  per  decade 
so  no  graph  is  available.)  The  statistical  significance  of  the  slope 
is  implied  by  the  b/Oj^  value  also  tabulated  for  each  curve.  The  var- 
iances of  range  we.e  also  averaged  ove**  all  SHOB  groups  pertinent  to 
each  range  and  the  resulting  standard  deviation  presented  graphically 
in  Figure  4.2-20.  Note  the  very  strong  linear  dependency  in  this 

composite  data  from  590  points.  Vhe  linearity  is  also  evident  in  the 

1/2 

sub-intervals  above  and  U^low  1 kft/kt  » the  approximate  scaled  range 

1/3 

to  10  psi  peak  overpressure  and  0.9  psi-sec/kt  ' overpressure  impulse. 
Thus  tlie  lognormal  distribution  is  implied  as  more  appropriate  than  the 
normal  distribution  for  both  the  peak  and  impulse. 
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Figure  4.2-11.  STANDARD  DEVIATION  OF  RANGE  VERSUS  R,ANGE  FOR 
0 to  3 HOB 
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Scaled  Range  (kft/kt*/^) 

Figure  4.2-12.  STANDARD  DEVIATION  OF  RANGE  VERSUS  RANGE  FOR 
5 to  11  ft/kt^''^  HOB 
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Figure  4.2-14.  STANDARD  DEVIATION  OF  RANGE  VERSUS  RANGE  FOR 
113  to  157  ft/kt’^^  HOB 
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Figure  4.2-15.  STANDARD  DEVIATION  OF  RANGE  VERSUS  RANGE  FOR 
182  to  205  ft/kt^/^  hob 
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Figure  4.2-19.  STANDARD  DEVIATION  OF  RANGE  VERSUS  RANGE  FOR 
1003  to  1249  HOB 
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Section  5 


ANALYSIS  FOR  BIAS 

The  potential  yield  adjustment  that  could  be  accomplished  by 
improving  the  random  error  variation  involves  the  question  of  how  much 
one  can  improve  the  yield  specification  if  the  damage  function  (e.g., 
a log  normal  distribution)  were  more  peaked.  The  potential  yield 
adjustment  from  removing  systematic  biases  could  result  in  changing 
the  yield  by  the  range-cubed  relationship  between  yield  and  range  to 
a specified  effect,  which  is  a stronger  relationship  between  yield 
and  range  variation  than  the  relationship  for  the  random  case.  Thus 
it  is  important  to  scrutinize  the  weapon  produced  environment  data 
for  possible  biases,  in  addition  to  assessing  the  effects  of  random 
errors. 

This  section  presents  the  results  of  analysis  of  the  nuclear 
test  blast  data  undertaken  to  ascertain  if  there  are  obvious  biases 
in  thesedata.  The  approach  for  this  analysis  was  to  identify  a set 
of  factors  that  could  potentially  cause  systematic  errors  and  to 
classify  the  data  according  to  these  factors.  An  analysis  of  variance 
was  then  conducted  on  the  data  to  indicate  whether  any  of  the  effects 
observed  were  “significant"  (in  the  statistical  sense),  or  could 
reasonably  have  arisen  simply  by  chance.  For  this  analysis,  three 
factors  were  chosen  to  discriminate  on,  namely,  instrument  type, 
yield,  and  wavefonr  type.  Three  levels  of  each  of  the  first  two  factors, 
and  two  levels  of  the  third  factor  were  defined: 

t Instrument  type:  SC  and  SRI  gauge  recordings;  NOL  instru- 

mentation recordings;  and  BRL  instrumentation  recordings. 

• Yield:  low  yield  (<3  Kt);  medium  yield  (3  Kt  to 
300  Kt);  and  high  yield  (>300  Kt). 

• Waveform  type:  near-ideal  (type  V);  precursor-associated 
(non  type  V). 

Only  tests  for  scaled  HOB  < 11  ft/kt^'^^  were  used. 


r 


I The  analysis  was  conducted  for  the  scaled  ranges  to  10,  30,  and 

i 100  psi  peak  overpressure  using  all  the  data  in  Reference  4.  A series 

j of  four  two-factor  cross-classified  experimental  designs  were  set  up 

5 

I and  analyzed.  Cross-classified  designs  discriminating  on  yield  and 

( instrumentation  type  were  constructed  for  scaled  range  to  10,  30,  and 

1 100  psi  data,  and  a cross-classified  design  discriminating  on  yield 

! and  waveform  type  was  constructed  for  scaled  range  to  30  psi  data. 

I Designs  discriminating  on  waveform  type  were  not  constructed  for  the 

i scaled  range  to  10  and  100  psi  data,  since  at  these  ranges  virtually 

! all  of  these  data  were  of  the  same  waveform  type  for  these  low  SHOB. 

I In  obtaining  the  ranges  to  the  nuclear  produced  environments  for 

i each  shot,an  attempt  was  made  to  interpolate  between  data  points  of 

{ a given  classification  so  that  the  resulting  estimated  range  to  each 

j peak  overpressure  was  based  on  a group  of  gauge  readings  rather  than 

I on  a single  reading.  In  the  case  where  only  one  instrument  type  was 

located  in  the  vicinity  of,  for  instance,  the  30  psi  range,  the  range 
was  estimated  from  the  "DASA"  curve  through  these  data.  In  the  case 
where  several  instrument  types  obtained  recordings  in  this  region,  an 
estimate  of  the  possible  bias  between  the  instrument  type  and  the 
1 fitted  line  was  obtained.  Thus,  if  an  instrument  type  showed  recordings 

consistently  higher  or  lower  than  the  fitted  curve  in  the  region,  this 
bias  was  accounted  for  in  estimating  the  range  to  30  psi  as  recorded 
by  that  instrument  type.  Figure  5.1  illustrates  tin’s  procedure.  For 
the  experimental  design  involving  waveform,  the  waveform  type  was 
taken  as  tabulated  in  Reference  4. 

I Table  5.1  shows  the  cross-classified  design  discriminating  on 

instrument  type  and  yield  for  range  to  30  psi.  The  data  found  to 
have  the  appropriate  yield/instrument  type  combinations  are  shown  in 

I 

each  cell,  along  with  the  estimated  range  to  30  psi  for  that  combination. 
Also  shown  are  the  cell  means  and  the  column  and  row  means  for  each 
factor  level.  The  analysis  of  variance,  showing  the  partitioning  of  the 
total  sum  of  squares  among  the  factors,  was  performed  on  these  data;  the 
i results  are  presented  in  Table  5.2.  Tables  5.4,  5.5,  5.7,  5.8,  5.10, 

; and  5,11  show  the  experimental  designs  and  analysis  of  variance  for  the 


for  the  10  and  100  psi  yield  versus  instrument  type  experiments,  and 
for  the  range  to  30  psi  yield  versus  waveform  type  experiment. 

To  ascertain  the  significance  of  the  variation  ascribed  to 
each  factor  and  the  interaction,  standard  F-tests  were  performed 
utilizing  the  mean  squares  from  the  analysis  of  variance  tables. 

These  are  presented  in  four  Tables,  5.3,  5.6,  5.9,  and  5.12.  As  can 
be  seen,  yield,  when  compared  to  the  residual  variatiOL,  appears  to 
make  a significant  contribution  to  the  total  variation  for  at  least 
3 of  the  experimental  designs.  In  fact,  for  the  range  to  30  psi  data 
in  the  design  discriminating  on  yield  and  instrumentation  type,  this 
contribution  can  be  described  as  “highly  significant"  since  a result 
like  the  one  observed  would  occur  by  chance  less  than  one  percent  of 
the  time  if,  in  fact,  there  is  no  true  contribution  due  to  yield.  For 
the  range  to  10  psi  data,  however,  the  result  appears  to  be  only 
marginally  significant.  In  fact,  if  it  were  not  for  the  other  results, 
all  of  which  indicate  yield  as  a factor  causing  a significant  amount 
of  variation  in  the  data,  the  results  for  the  range  to  10  psi  data 
would  have  been  dismissed  as  a "real"  effect,  since  results  of  this 
type  could  be  observed  a sufficiently  large  number  of  times  by  chance 
to  discount  the  results  as  significant. 

To  ascertain  the  possible  effect  on  the  change  in  yield  suggested 
by  the  yield  scaling  bias,  consider  the  case  for  the  range  to  30  psi 
data: 

• The  mean  of  the  medium  and  large  yields  (computed  from 

Table  5.1)  is  575.72  ft/ktl/3  for  30  psi.  This  represents 
a systematic  bias  between  small  and  1a>'ge  yields  that 
corresponds  to  an  error  factor  of  f^  = 1.08,  which  in 
turn  would  correspond  to  a yield  mis-specification  of 
about  27  percent  if  the  scaled  results  of  large  yields 
were  used  to  represent  small  yields.  (The  same  argument 
applies  if  the  scaled  results  of  small  yields  are  to  be 
used  to  represent  large  yields.) 

e The  scaled  range  to  30  psi  from  DASA  2506  is  610  (572 
based  on  using  a multiplication  factor  of  only  1.65  for 
the  yield  "doubling"  expected  at  the  surface).  If  this 


is  used  to  represent  1c#w  yield  environments,  a mean 
bias  corresponding  to  an  error  factor  of  fjj=1.15  (1.07 
for  1.65  "doubling")  is  indicated.  This  corresponds  to 
a mis-specification  of  yield  by  approximately  50  percent 
(24  for  1.65  "doubling"'). 

The  above  suggests  that  there  might  be  a relationship  between 
yield  and  the  scaled  range  to  30  psi  for  these  surface  tests.  In 
other  words,  cube-root  scaling  may  not  be  appropriate  over  the  entire 
yield  spectrum.  Figure  5.2  shows  a plot  of  yield  versus  scaled  range 
to  30  psi;  the  line  through  the  points  was  fit  to  all  15.  Visually, 
there  appears  to  be  a linear  relationship  between  the  log  of  yield 
and  scaled  range,  i.e.,  a bias  that  exists  even  after  the  initial 
scaling  for  surface  shots.  Another  observation  is  that  Sugar  (1.2  kt) 
and  Koon  (150  kt)  appear  to  be  outliers  on  this  plot.  Sugar  was  an 
early  test  conducted  with  early  instrumentation  (but  it  was  found  that 
Instrumentation  did  not  contribute  to  the  variation),  and  Koon  was  a 
test  in  which  significant  instrumentation  "failures"  were  observed 
(although  DASA  1200  reports  that  the  recordings  presented  are  good 
data).  No  immediate  explanations  are  available  as  to  why  Sugar  and 
Koon  should  be  outliers. 

For  the  present,  assume  that  there  is  a linear  relationship 
between  the  log  of  yield  and  scaled  range  to  30  psi,  and  assume  as  a 
tentative  hypothesis  that  this  phenomenon  is  due  to  a variation  in 
the  fraction  of  total  energy  that  goes  into  the  nuclear  produced 
blast  environment  that  scales  with  yield.  Then  the  true  scaling 
relationship  between  yield  and  range  would  be  the  form: 


where 

= True  range  to  30  psi  for  arbitrary  yield  W 
Rq  = Range  to  30  psi  for  1 Kt 
Wq  = 1 Kt 
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Table  5.2 


Sum  of  squares  of  Factor  B 


F-TESTS  FOR  SIGNIFICANCE 
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Table  5.4 

YIELD  VS.  INSTRWtEMTATION  TYPE  EXPERIKENTAL  DESIGN 
SCALED  RANGE  TO  ID  PSl 


VIELD  VS.  INSTRUMENTATION  TYPE  ANALYSIS  OF  VARIANCE 
SCALED  RANGE  TO  10  PSI 


Table  5.6 

YIELD  VS.  INSTRUMENTATION  TYPE  F-TESTS 


INSTRUMENrATION  Jm  EXPERIMENTAL  DESIGN 


YIELD  VS.  INSTRUMENTATION  TYPE  F-TESTS 
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i 

i 

i ^0  “ multiplicative  factor  for  1 Kt  (proportional 

i to  energy  coupling  into  air  blast) 

g(W)  = Yield  multiplicative  factor  for  yield  W. 

From  Figure  5.1,  the  ratio  of  to  should  be  "corrected"  by 
the  factor  S^R^,  where  is  the  scaled  range  to  30  psi  calculated 

from  the  (W)  1/3  relationship  (i.e.,  the-"old"  scaled  range).  Thus: 


but 

“ Ro  k log  W (from  Figure  5.1)  . 

Therefore , 


where 

k'  - k/Ro  . 

Therefore,  the  yield  multiplicative  factor  related  to  energy  coupling, 

g(W).  is: 


g(W)  « l/(l+k’  log  W)^  . 


To  evaluate  k*,  use  the  expression  for  S^; 

= Ro  I09  w 


Using  Figure  5.2  and  20  Kt; 

. 566  - 552 

“ ■ 1.301 

= 10.76 

k'=  k/Rg 

= .0195 

The  semi -empirical  expression.  Equation  (5-1),  can  be  used 
to  scale  low  HOB  tests  and  should  be  adequate  between  10  t and  10  Ht. 
This  correction  is  applicable  only  for  30  psi  from  surface  bursts. 

The  potential  biases  for  scaled  ranges  to  10  and  100  psi  are 
not  as  pronounced  as  the  bias  for  the  range  to  30  psi;  in  fact,  there 
appears  to  be  no  readily  discernible  bias  characterized  by  the  low 
yield  test  results  being  significantly  different  from  the  medium  and 
high  yield  test  results  for  the  scaled  range  to  10  psi  data  (the 
analysis  of  variance  indicated  a less  significant  result  concerning 
the  amount  of  variation  in  the  data  correlating  to  the  factor  yield 
for  this  nuclearly  produced  environment  also).  The  implication  is 
that  the  effects  of  blast  prooagation  override  any  second-order 
effects  that  correlate  to  the  factor  yield  at  the  scaled  range  to 
10  psi. 

For  comparison.  Figures  5.3  and  5.4  show  plots  of  yield 
versus  scaled  range  to  10  and  100  psi,  respectively. 
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Appendix  A 

THE  LOGNORMAL  DISTRIBUTION 


A random  variable  x is  said  to  be  lognormally  distributed 

2 

if  its  logarithm  is  normally  distributed  (s  = variance  of  In  x) 


dp  _ J. 


exp 


n(x)  - ln(Xg 


2s 


I 


ln^(x/Xp) 


= L (x:  x^.s)  and  o<x,Xjj,s, 


(Al) 


The  lognormal  distribution  is  a natural  distribution  in  that  it  is 
appropriate  for  any  random  variable  which  is  the  product  of  random 
variables.  In  fact,  under  certain  not  very  restrictive  conditions, 
a Central  Limit  Theorem  can  be  proved,  stating  that  as  N goes  to 
infinity,  the  product  of  N independent  random  variables,  no  matter 
how  they  are  individually  distributed,  is  lognormally  distributed 
(References  1,  2}.  Furthermore,  the  random  variable  defined  as  the 
sum  of  lognormal  variates  is  in  many  cases  accurately  represented 
by  a lognormal  distribution  instead  of  the  normal  distribution  that 
might  be  expected  on  the  basis  of  the  Central  Limit  Theorem  (Reference  3). 
Therefore,  if  a random  variable  can  be  written  as  the  sum  of  products 
of  random  variables,  in  many  cases  the  appropriate  distribution  is 
the  lognormal.  For  very  small  variance,  the  region  around  the  mode  of 
the  lognormal  can  be  approximated  by  a normal  distribution  but  the 
tolerance  limits  cannot  be  as  easily  approximated. 

Many  properties  of  the  lognormal  distribution  can  be  explicitly 
derived,  and  will  be  listed  here  for  use  in  Section  3. 

The  cumulative  distribution  function  (CDF)  is  given  in  terms 
of  the  normal  CDF: 

dx  = F I ^ 1 where  F(z)  = [l  + erf  (z//T)|  (A2) 
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The  median  of  the  distribution  is  thus  x.  from  z=0.  Most  statistics  of 
0 

interest  can  be  derived  from  the  moments  about  the  origin: 


al  <«) 

= X."  exp 


"0  ' ' 2 

so  that  the  mean  and  variance  of  x are  given  by 

-2 


(A3) 


X = 


7'=  “o  (t) 


= x2  - x^  = x^  (exp  (s^)  - 1) 


and  s^,  the  logarithmic  variance,  is  given  by 


(A4) 


s^  “ In  (1  + ( o'/x)^). 


(A5) 


i 


# 

The  lognormal  has  a very  simple  expression  for  the  moments  of  x about  the 
origin;  the  expression  for  moments  about  the  mean  is  more  complex. 

This  is  just  the  reverse  of  the  case  for  the  normal  distribution. 

The  mode  of  the  distribution,  where  dp/dx  is  maximum,  is 

given  by 

’‘mode  = ’‘o 

The  lognormal  distribution  is  always  skewed,  with 

W<V^  (''7) 

This  is  illustrated  in  Figure  A-1  for  unit  median  and  logarithmic 
variance. 

The  percentiles  or  m-tiles  of  the  lognormal  distribution, 
x,^,  are  defined  as  usual 


A- 2 


1 


so  that  Xq  = XgQ.  A synrietric  p confidence  interval  is  conveniently 
defined  around  x.  in  terms  of  a confidence  level  error  factor  f_ 

0 p 

such  that 


x>.  = f„x^ 

and 

m p 0 

’‘lOO-m  = 

0 p 

where 

(A9) 

m 

I(30  = (p+l)/2 

. 

Then  it  can  be  seen  that 


(AlO) 


where  F is  given  by  Equation  (A2).  This  allows  the  parameters  x^.s 
to  be  determined  entirely  from  the  error  bounds  on  the  confidence 
level  (CL),  i.e.,  x^jj,  Xjqq_j|j.  Typically  this  is  done  for  a 90  percent 
CL:  p=0.9,  m=95.  For  this  choice  the  notation  in  this  report  omits 
the  subscript  p=0,9  from  fpt 

f^  = Xgg/Xg  (905iCL) 

*0  V ^95*5 


Two  important  theorems  about  the  lognormal  distribution 
are  used  in  this  report.  The  first  is:  _ 

Theorem  1:  If  x is  a lognormal  variate  with  parameters 
X-,s,  then  z = ax^  is  a lognormal  variate  with  parameters 

O u 

axg  , bs. 

Theorem  1 is  proved  by  direct  substitution. 
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Theorem  2;  If  x,  y are  lognormal  variates  with  and 

^o’®y’  ^ ~ * lognormal  variate  with  Zg^x^^’^y^'^ , 


2 2 2 2 2 
=2  = P 'x-^^  S- 


Theorem  2 is  a standard  theorem  proved  many  places,  e.g.. 

Reference  4,  p.  278  ff.  However,  the  proof  is  illustrative  of  some 
of  the  manipulations  required  in  this  study,  so  will  be  briefly 
outlined  here.  From  Theorem  1,  xP  and  y*^  are  lognormal  variates, 
so  Theorem  2 is  similar  to  the  theorem  that  the  sum  of  two 


normal  variates  is  a normal  variate.  For  this  purpose  define 
normal  variates  u,v  with  means  Uq,v^  and  variances  o^,  o^,  and 
a random  variable  w,  with  mean  w^,  a^.  Then  the  probability  that 
the  sum  of  u and  v is  less  than  some  particular  value  w is  given  by 
the  probability  of  occurrence  of  a particular  u,  together  with  the 
probability  that  v is  less  than  or  equal  to  w-u  for  all  u less  than  w. 
This  relationship  may  be  expressed 


1 

j 

i 

/.w 

-J 

/•W 

f ^ 

f du 

^-u 

1 ^ dv  du 

(A12) 

1 or 

•f 

\ 

i ; 

j ; 

Pw(w)  =j 

' du 

p^(w-u)  du 

(A. 3) 

It  is  convenient  to  define  the  Laplace  transform 

where^is  the  Laplace  operator  and  not  the  PDF  of  Equation  (Al). 
Then  using  the  common  properties  of  the  Laplace  transform,  the 
convolution  in  Equation  (A13)  becomes 

i!{Pw.s)  =i!(^.s)i!(p^,s)  (A15) 
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and  the  relation 


s j;^(p„.s)  =£  (^.  s)  (A16) 

gives  the  result 

£(^.  s)  s)4;(f . s)  (A17) 

It  is  easy  to  see  that  for  the  normal  distributions.  Equation  (A14/ 
gives 


= exp  (-  u^s  +<r^^V/2).  (A18) 

As  a result  the  unknown  dp/dw  has  the  Laplace  transform 

^ ^ Ov^)s/2),  (A19) 

which  is  the  Laplace  transform  of  a normal  distribution  with  mean 

2 2 2 

w„  = u„  + v„  and  variance  o,,  = o + o , For  the  present  case 
0 0 0 w u V 

u=tn(xP),  Oy^PSj^  . with  similar  relations  for  v.  This  substitution 
proves  the  theorem.  Theorem  2 can  obviously  be  expanded  to  include 
random  variates  that  are  the  products  of  any  number  of  lognormal 
variates. 
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